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Abstract 

Due to relativistic effects, the observed emission from relativistic jets Is 
quite different from the rest frame emission. We discuss systematic differences 
between the observed and intrinsic Intensities of sources In which jet phenomena are 
occurring. Assuming that jets have a power law luminosity function of slope B, we 
calculate the observed luminosity distribution as a function of the velocity of the 
jet, the spectral Index of the rest frame emission, and the range of angles of the 
jets relative to our line of sight. The result Is wel 1 -approximated by two power 
laws, the higher luminosity end having the original power law Index B and the lower 

luminosity end having a flattened exponent Independent of B and only slightly great- 

er than 1. We then Investigate a model consisting of beamed emission from a jet and 
unbeamed emission from a stationary central component. The luminosity functions for 
these two-component sources are calculated for two ranges of angles. For sources In 

which beaming is Important, the luminosity function Is much flatter. Because of 

this, the relative numbers of "beamed" and "unbeamed" sources detected on the sky 
depend strongly on the luminosity at which the comparison Is made. 

mso Johns Hopkins University. 

^Now at Institute of Astronomy, Cambridge, England. 
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I. Introduction 


OF. PXS^M QUAUTf 

Relativistic jets emanating from the cores of compact extragalactic objects 
have been Inferred from high resolution radio maps, superluminal expansion, and 
brightness temperature limits. Roughly 70% of the VLA- mapped radio sources have 
jets, some of which may be relativistic (Bridle 1982). If these relativistic jet 
phenomena are common, the relativistic motion of the emitting plasma has a tremen- 
dous effect on the observed emission and the consequences for observational work are 
striking. Specifically, the observed luminosity of a relativistically moving source 
may be. considerably brighter or considerably dimmer than the luminosity emitted In 
the rest frame, depending on the direction of the motion with respect to the obser- 
ver. Taking this effect Into account Is Important In the derivation of luminosity 
functions, In the calculation of energy conversion efficiencies, and In any discus- 
sloe of the physics of sources with relativistic jets. 

This paper Investigates the difference between observed and intrinsic lumino- 
sity functions for sources In which relativistic jets are important. A two-compo- 
nent model, consisting of a > st plus a stationary component, Is considered. When 
the jet points within a certain angle, 6^,, of the 1 Ine-of-sight, Its relativistical- 
ly enhanced emission will dominate; outside of this angle, the jet is decreasi ( 
Important. We calculate luminosity functions for these two-component sources for 
two ranges of jet angles; O®<0<0^, and 0j,<e<9O°. The results are applicable to any 
frequency band, assuming only a power law continuum with energy index a. (For a 
more complicated spectral shape, the power law approximation will still be valid for 
a sufficiently narrow bandwidth. The implications for sources with complex spectra 
are discussed in §111.) A forthcoming paper, Urry and Mushotzky (1984), discusses 
the relevance of the quite general result presented here to BL Lacertae objects, in 
which the jet may be the distinguishing characteristic. In §11 of the present work 
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we describe our assumptions and outline the calculation. In §111, results for 
representative parameters are shown and the Important features are enumerated. In 
§IV, we summarize our findings and discuss the Implications. 


II. Luminosity Enhancement 


ORIGINAL OTi’ 
OF POOR QUAUW 


The observed luminosity, L, of a relativistic jet Is related to Its emitted 
luminosity, ^ via 

L = 6P,jC, (1) 


where 6, the kinematic Doppler factor for the jet, 1$ defined by 6= (y(1-0 cos0))"1, 
3=veloc1ty in units of the speed of light, 1s the Lorentz factor, and 0 

Is the angle between the velocity vector and the line-of-sight. 

The exponent p depends on assumptions about the spectrum of the jet emission 
(Ryle and Longair 1967), the structure of the jet (Scheuer and Readhead 1979), and 
the frequencies at which the comparisons are made. For the spectrum we assume a 
power law of the form F^~v"“. Then If L refers to observed monochromatic luminosity 
(ergs s“i at frequency v compared to .Jt{v), p=3+a. Equivalently, L(vg)= 6 ^ 

^(vg), where the frequency of observation is 6 times the emitted frequency Vg. 

Two powers of the enhancement are due to aberration: the emission is beamed 1n the 

forward direction due to the relativistic motion of the jet. Another power comes 
from contraction of the time interval by a factor 6 so that the number of photons 
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per unit time 1s enhanced by The remaining power of a present In the comparison 
at a fixed frequency Is due to the actual blue-shifting of the photons. (The pho- 
tons observed at frequency v were emitted at frequency so that there were a 
factor 6°^ more than were emitted at frequency v.) If Instead L refers to integrated 
luminosity (ergs s”M over corresponding bands, then P'^4, since the blue-shifting of 
the photons changes the observed bandwidth by a factor 6. Scheuer and Readhead 
(1979) noted that if discrete components are ejected, 1n pairs, the observed life- 
times of approaching and receding components differ, with more of the receding blobs 
being visible at any one time. Under these conditions, and with reference to flux 
densities or monochromatic luminosities at a fixed frequency, p=2+a. In practice, 
the receding jets are so underlumlnous due to relativistic effects that In the 
context of a two-component model (jet plus stationary source), they may be Ignored. 
Calculations for 90°<e<180° show that even if the beamed and unbeamed components are 
comparable in Intrinsic luminosity, the receding beamed component will be less than 
one percent of the total luminosity for likely values of y. 

The larger the value of p, the greater the effect on the observed luminosity. 
Thus, sources with intrinsically steep spectrum emission are more affected than 
those with flat spectra. If the source spectrum is more complex than a single power 
law, its shape remains unchanged but (on a plot of flux density versus frequency, 
and with 6>1) it Is simply shifted up by a factor of 6^ and to the right by a factor 
of 6. Observations in a given band therefore depend crucially on what is emitted at 
frequencies a factor of 6 lower, so that complex spectra can give rise to serious 
complications. A case In point Involves X-ray observations of BL Lac objects, which 
may be dominated by relativistically beamed emission (Wei strop 1981; Urry and 

Mushotzky 1982; Worrall and Bruhweiler 1982; Urry 1982a, b; Worrall et al. 

1984). The few BL Lac X-ray spectra that have reasonable statistics suggest there 
are two components to the spectrum, a steep one at low energies and a flat one at 
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higher energies (Mushotzky et a1* 1978; Riegler, Agrawal, and Mushotzky 1979; 

Worrall aU 1981; Urry and Mushotzky 1982; Urry aK 1982b). The location of 
this steep- to-fl at break seems to occur anywhere from the extreme ultraviolet to the 
medium X-ray (~10 keV) range. Thus, if the X-rays originate in a relativistic beam, 
the degree of enhancement in the X-ray band can be determined only with sufficient 
spectral information. In this paper we have Ignored the spectral complications and 
simply assumed the emission is described by a simple power law. 

We begin with a large number of sources, each having a pair of oppositely 
directed jets so that each source has an approaching jet. Our basic assumption is 
that these jets are randomly oriented on the sky, i.e. the probability distribution 
in e is uniform («<d{cose)). For simplicity we assume that all jets have the same 
fixed Lorentz factor. While not exact, this is a good approximation if the range in 
Y is small. The results presented below used a y of 5, as suggested by the results 
of Orr and Browne (1982), and likely values are in the range 2 to 10 (Blandford, 
McKee, and Rees 1977, and references therein; Cohen _et_£l_. 1977; Marscher and Brode- 
rick 1981; Sherwood e^^. 1983; Unwin al_. 1983; Readhead £t 1983; Madejski 
and Schwartz 1983; Henrikson, Marshall, and Mushotzky 1983). For 2<y<10 the shape of 
the observed luminosity function is qualitatively the same: the normalization is 
proportional to y”^ and the luminosity at which the two power laws join is propor- 
tional to Y*^ (see Eqn. 7 below). Thus, introducing a distribution in y would tend 
to blur the sharp two-power law form of the resulting observed luminosity function, 
smoothing the break between the two power laws and spreading it out over a range of 
luminosities. The true distribution of Lorentz factors in extragalactic jets is 
unknown, but our results for the extreme values of y will describe the envelope of a 
more realistic calculation. 

Since the Doppler factor, 6, is a function of e and y only and the latter is 
fixed, the probability distribution in 6 is easily derived from P(e): P(6)= 
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(0y 62)’"1, jg twice the probability calculated by Schwartz^ Madejski, and Ku 
(1981), corresponding to the picture that jets come in pairs. Given a particular 
emitted luminosity we then know the probability of observing luminosity 
L« 6^ dl* That is, 

P(LW) = P(5)M = ^ . (2) 

dL 0YP 

For any form of intrinsic differential luminosity function the observed 

differential luminosity ^unction will be given by 


4,o(L) = / di4,e(i) P(LlX) . 


(3) 


iSaVw I JCvwu vuf W 


simple luminosity function? the form 





or ^>^ 2 * 


(4) 


The sharp cutoffs are almost certainly unrealistic but they conform to our practice 
of using the simplest picture that is still consistent with observation, in this 
case that luminosity functions must turn over at low luminosities (Piccinotti et al . 
1982; Elvis, Sol tan, and Keel 1983; Gaston 1983). 

Combining equations (2) and (4) with equation (3), we can perform the integral 
over to obtain the transformed luminosity function, It will be useful to 

consider the luminosity function for various fixed ranges of 0. Let and 
represent the limits of the range under consideration. Then 6 has minimum and 
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maximum values corresponding to and respectively. Clearly 6,^gj^a (y( 1- 
$cos0,^j^))"i and (Y(l-0cos0 j„ax^ limits 0*^ and 90^", «max“2v and 

Thus, only part of the L-X plane Is accessible, as Illustrated In Figure 
1. The limits of the Integral depend on the value of Lr for which there are three 
regions. These are indicated In Fig. 1 by the dashed lines at 5|T,in^<>Ci» >-3“ 
^m1n^^2» ^4"* ^i» ^m)C ^max^'^a’ result for the observed 

luminosity function has two forms, depending on whether is less than or greater 
than L^. With the abbreviation C« B-l/p-1, the results are: 


PACS IS 

QUAUTt 



^^Uln 

*-m1 

L3<L<L^ (5) 



and when L3>1^, 
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0 , 


L<L 


min 


,-C -(p+D/p pC -8 

L “ °m1n ^ * 




♦od) = ^ < l'® (Cx - 

pYPL» 


L^<L<Lg (6) 


pC -B .-C -(pH)/p 

'^max ^ <2 ^ 


L3<L<L^ax 


0 , 


^^^-max* 


In each Hne, the first term dominates the second, except at the limits. Neglecting 
the second terms, both cases reduce to the same approximation: 


/ 


l<L 


min 


c|>o(L) - Lmin<L<^ 


(7) 


6^^ l“® 
°max ^ » 


V 


0 , 




*-^*-msx 


Thus the beamed object luminosity function is a broken power law. At the low lumi- 
nosity end the index is (p+l)/p, which is between 1 and 1.5 for reasonable values of 
p. All of the objects in this range have the lowest intrinsic luminosity and 
are Doppler boosted by the full range of Doppler factors, so that the observed 
luminosity function <))q(L) simply follows the probability distribution of Eqn..(2). 
Toward the high lumiryosity end there is a break at L= 1^^= ^mai< ^ ^ (steeper) 
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power law with the original Index B, This occurs because the large number of sour- 
ces at the low luminosity cutoff cannot be boosted any higher than L^, so the number 
of sources observed above that luminosity falls off exactly as the number of sources 
existing above although the normalization Is different. If the lower cutoff Is 
not sharp tut a gradual rollover, the break L^^ will be similarly gradual. The 
strong dependence of the result on the low luminosity cutoff Is explicit In Eqn. 

(7); the upper cutoff plays no role at all, except for determining Lj^g^. In prin- 
ciple, It Is possible that the moderate luminosity beamed sources from the flat part 
of the observfid luminosity function may arise from a population of sources Intrinsi- 
cally too faint to have been seen when unbeamed. 


III. Results 

A representative calculation Is shown In the two panels of Figure 2. We as- 
sumed that the jet was a continuous stream of outward-moving particles (Konigl 1981, 
Marscher 1980), and we considered luminosity comparisons made at a fixed frequency, 

V that p*3+a. The functions plotted In Fig. 2 used the value p=4, corresponding to 
a spectral Index a=l. Radio spectra are often flatter (Scheuer and Williams 1968, 
Kellerman and Paullny-Toth 1969, Condon et_^. 1978); optical, ultraviolet, and X- 
ray spectra have a wide range of spectral Indices {Richstone and Schmidt 1980, 
Mushotzky et al. 1980, Worrall 1981, Urry et al. 1982a, Bregman £t 1982) 

but a~l Is a good representative value. We assumed all of the luminosity was In the 
jet, and the value used for the Lorentz factor was v=5, as suggested by the results 
of Orr and Browne (1982). Values for e were allowed to range between 0° and 90® 


9 


ORIGINAI FAGS Si 
OF POOR QUALITY 

stnce we are neglecting radiation from jets which point away. We defined a critical 
angle such that 6 (y, 6^;H. With is ^35?3, 

In Figure 2 we have plotted log CL.})q(L) 3 versus log L. The solid line in the 
top panel represents the intrinsic luminosity function used (Eqn. 4), with slope 
B»2.75 and cutoffs and (The units are arbitrary but a tic mark on 

either axis represents a decade.) The dash-dot line in Fig. 2a shrws the observed 
luminosity function as it has been altered by the effects of relativistic beaming. 
All angles are included, 0®<e<90®, corresponding to Eqn. 7 with 6|^j^p6(90®) and 

2b the obdects have been separated according to their orienta- 
tion angle, with the dotted curve representing the luminosity function for jets with 
angles 35?3<6<90® and the dashed curve representing jets with angles 0®<$<35?3. As 
expected, small orientation angles correspond to high observed luminosities. 

As discussed above, we do not expect all of the luminosity of a source to arise 
in a jet. Instead, let the total luminosity of a source Lj be composed of an un"* 
beamed part and the jet luminosity Lj. The simplest assumption we can make is 
that the intrinsic luminosity of the jet is some fixed fraction f of the unbeamed 
luminosity, i.e., Then 

Lt = otu + Lj 

= (1 + ffiP) . (8) 

The probability P(L|3C), where L=^T refers to the total luminosity and refers 

to the unbeamed luminosity, is again derived from P(6) according to Eqn. (2), with 
the resul t 

P(LtX)=_l_ i"’' (^- . (9) 

The Integral in Eqn. (3j is now done numerically and the results are shown in Figure 
3. Once again the luminosity functions have been separated according to orientation 
angle. In this case e^, was defined as the angle at which the unbeamed luminosity 
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was equal to the beamed luminosity of the jut* That Is, f<sP{0jj) Is I, or 0^« 
arcos(0'"^ (1 - Y"if^/P)). Thus, the critical angle changes as the fraction f Is 
varied. 

Figure 3 Illustrates the results for this two-component model, using the same 
Intrinsic luminosity function as In Figure For each value of f the sources are 
divided Into two classes; the jet-dominated sources with O®<e<0g and the non-jet- 
domlnated objects with 0g<0<9O®. The dashed curves In Fig. 3 represent the lumino- 
sity functions for jet-dominated sources. The values of f shown are 0.001, 0,01, 

0.1 and 1.0, for which the corresponding critical angles are 10?2, 17?0, 25?1, and 
35?3, respectively. The solid curves represent the luminosity functions for the 
non-jet-dominated sources. In these sources the jets point outside so that 
relativistic effects actually diminish the observed jet Intensities. The four 
curves are nearly Identical, lying one on top of the other, and are effectively the 
same as the unbeamed, Intrinsic luminosity function assumed. This Is the case both 
because the jet luminosity 1s insignificant when fsP Is less than unity, and because 
only a small number of objects are jet-dominated and are subtracted from 

the initial population. 


IV. Discussion and Conclusions 

We have described the impact of relativistic enhancement on measurements of the 
luminosity function for objects with relativistic jets. We divide such sources into 
two classes; those which are jet-dominated and those In which the jet contributes 
negligibly to the total luminosity of the source. Let the latter objects be called 
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parents, In reference to the parent po9fl?69R SH^^Pfves rise to the class of jet- 
dominated objects. The main results of our study are illustrated in Figure 3. 

First, the jet-dominated objects tend to be more luminous than the parents, although 
there are fewer of them overall. 

Second, there are luminosity ranges within which the numbers of jet-dominated 
objects are comparable to or even greater than the numbers of parents. For example, 
in Figure 3 the numbers of jet-dominated objects and parent objects for f«l,0, 0.1, 
and 0.01 are rougly equal at La8.6(ip L«13cfj* and L»20oifj, respectively; for 
fsQ.OOl, there are more parent objects at all luminosities below This is in 
marked contrast to the naive expectation that because Q<6j. only Infrequently, we 
should see far fewer beamed objects than parents irrespective of luminosity (or f 
for that matter). 

Third, over a large range in luminosity, the jet-dominated objects will have a 
flat lumir»i>;ity function, generally much flatter than is typical for suggested 
parent K.pufations. The flat slope results from the many objects at the lowest 
luminosity being boosted by a range of Doppler factors. Thus we are making a speci- 
fic prediction: if a given class of object consists of sources with observed emis- 

sion primarily from relativistic jets pointed at the Earth, that class should have a 
flatter luminosity function than the class of similar sources with jets pointed in 
other directions. This may explain the apparently flat slope of the X-ray luminosi- 
ty function for BL Lac objects (Schwartz and Ku 1983) compared to the slope for any 
Other class of extragalactic object (Piccinotti £t 1982). Because of the flat- 
tening of the luminosity function for jet-dominated objects, the ratio of the volume 
density of the jet-dominated objects to the volume density of the parents varies 
with observed luminosity. Estimates of the required volume density of the parent 
population giving rise to BL Lacertae objects have not taken this effect Into ac- 
count (Schwartz and Ku 1983, Browne 1983), 
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Jets jre present In extragalactic objects* If the jets are rslat1vfst1c> the 
proper treatment of the luminosity enhancement effect Is Important. If 3C273, for 
example, exhibits superlumlnal expansion because of a directed jet, (Cohen et al. 
1971, 1977, 1979; Sherwood et al. 1983) It may be much less luminous Intrinsically 
than observed, at least In the radio, depending on the value of f for the source. 

If the jet contribution to the total luminosity Is significant, one should not 
expect to see a few dozen sources as bright as 3C273 1n the same co-movIng volume as 
their jets will likely be pointing away. Finally, If this beaming picture Is cor- 
rect, the highest luminosity objects observed should be those which are dominated by 
a relativistic jet. These should tend to have strong nonthermal (jet-dominated) 
continua, and to be highly variable, highly polarized, and possibly superlumlnal . 
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Figure Captions 
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Figure 1, Observed luminosity L versus emittel luminosity li. 

^2 cutoffs for the rest frame luminosity function and and 

are the resulting cutoffs In the observed luminosity function. The dashed lines 
Indicate the ranges for L, which are delimited by the values !-niin"%!i ^ 3 "'^ml 8 
and l-j|,gj^“ 6 f,^g^ c^se L 0 <L^ Is Illustrated. In the two-compo- 

nent model, where the Intrinsic jet luminosity Is’filf the unbeamed luminosity Is 
the expression (l+ffi^P) must be substituted for 6 jP, with 1 =m 1 n or max. 


Figure 2. Luminosity Functions for Jet Component. 

(a) The solid line represents where <|)g(3C), the Intrinsic differential lumi- 

nosity function assumed, Is a power law given by Eqn. (4) with index B-2.75 and with 
luminosity cutoffs and 0 ^ 2 = 10 '*. The units are arbitrary. The dash-dot line 
Is the observed luminosity function due to relativistic beaming with y= 5 and p=4, 
for sources at all angles 0°<e<90®. 

(b) The observed luminosity functions for beamed sources are separated into two 
groups according to jet orientation angle. Dotted line; jets with angles 
35°3<0<9O°. Dashed line: jets with angles 0°<6<35^3. The sum of these two lines 
is given by the dash-dot line in panel (a). 


Figure 3. Total Luminosity Functions for Two-Component Model with Beamed and 
Unbeamed Emission. 

The observed luminosity functions for sources with total luminosity comprised of the 
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unbeamed luminosity and the beamed (jet) luminosity These sources have 

been divided Into two classes: those dominated by unbeamed emission as defined by 
0^,<0<9O® (solid lines), and those dominated by beamed emission from a jet with 
O®< 0 < 0 ^. (dashed lines). Four values of f were used, as Indicated in the figure. 

The curves for non-jet-domlnated sources are virtually unaffected by the value of f, 
and depart from a single line only at the endpoints. 

OF POOR 
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